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Smoothing rough edges 


Dan O’Brien spent last month figuring out how to make 
cartoon characters walk. This month he enlists Bezier’s 
help for a new challenge: to make them walk smoothly. 


Last month, Low Level wentinto some 
detail to explain a program which cal- 
culated the leg joint angles of a walking 
human. Fun to watch, but not the most 
universally applicable of algorithms, 
unfortunately. This month, as an ad- 
dendum to WALK.BAS, we’ll concen- 
trate on a set of procedures that are still 
fun to watch but have a more wide- 
spread application: curve drawing func- 
tions. These can generate smooth, 
natural looking outlines to graphics, 
text characters and — in the case of 
WALK.BAS — graphs of functions. 
We left WALK.BAS with one niggle. 
While the maximum and minimum 
values of each joint angle was correct, 
the movement of the joint itself be- 
tween these two numbers was some- 
what jerky. In the manner of some robot 
arms (and all robotic mime artists), the 
leg joints move at an even speed until 
they reached a prescribed angle, then 
stop dead. Immediately afterwards, the 
joint begins the move to the next set 
angle, at a rate bearing no resemblance 
to its speed ofarrival. Plotted as acurve, 


the function of time to angle would be 
a set of spikes, with the points of the 
spike marking our set maxima and 
minima. 

Animals, humans, and indeed most 
of nature does not work like this. Gen- 
erally, changes in the natural world are 
described by curves. So, a human knee 
tends to slow its speed of turn as it 
reaches a maximum. But, oh, how can 
a poor and humble programmer ape 
nature’s beautiful and languid poise? 

Easy. He or she steals. In this case, 
the best people to steal from are the 
mathematicians. Many of them are very 
weedy, and are quite happy to lend a 
few concepts. If they seem too obtuse to 
nick from, try geometricians. Ata push, 
even engineers will do. 

The first curve drawing algorithm 
designed for programming use was pil- 
fered wholesale from an engineering 
toolbox by Pierre Bezier, a programmer 
on one of the first computer aided de- 
sign projects at Renault in the mid- 
sixties. The tool he stole was the spline. 
This long strip of flexible metal could 
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CHANDI SOFTWARE 
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HARROW, MIDDX. HA1 2NB 
SOFTWARE NOW AVAILABLE AT DISCOUNT PRICES 


Tel: 081907 6497 
Fax: 081 907 3904 


Amstrad DMP 2000/3000 
Amstrad DMP 4000 

Citizen 120D/LSP10/Swift 24/9 
Citizen Swift 24 Colour 
Commodore MPS1220/1230 
Epson LQ400/500/800/850 

Epson FX/MX/RX80/FX/LX800 
Epson FX/MX/RX100/FX/MX1000 
Epson LX80/86/90 

NEC Pinwriter P2200 

Panasonic KXP1123/1124/1140 
Panasonic KXP1080/1180/90/1592 
Star LC10/20/100 Black 

Star LC10/20/100 Colour 

Star LC200 Black 

Star LC200 Colour 

Star LC24-10/200 Black 

Star LC24-10/200 Colour 

Canon BJ-10 Inkjet Cartridge 
HP Deskjet Cartridge (Double Cap) 
Canon BJ-10 Refill (Twin Pack) 
HP Deskjet Refill (Twin Pack) 


Ring for ribbons and inkjet cartridges/refills not listed 
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mouse. Faders page, select sounds by name, 
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Postscript printers & Postscript fonts (for 
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Demo version - £5 deposit Turtle Beach Multisound (Itd offer) £499 
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SCIENTIFIC SHAREWARE 


We have the largest, range of technical, scientific 
and hard to find Public Domain and Shareware 
for the IBM PC in Europe, including:- 


MATHS & STATS, MEDICAL, ELECTRONICS, HAM 
RADIO, ENGINEERING, ESOTERIC, REFERENCE 
DATA, LANGUAGES & PROGRAMMING, SOURCE 

CODE, GRAPHICS, WINDOWS, BUSINESS and lots more. 


PDSL is the longest established library in the UK, 
Since 1982 we have been acknowledged as wings 
the best most reliable library service in the UK. 


If i too would like to find out about the true wealth 
of PD and Shareware available, phone, write or FAX 
today for a free prinied catalogue, 

(We also stock many registered versions.) 


fo] The Association 
O of 


Our unique system for 


Wi Writing 


Starts you in business immediately hadaue. Macrhees hte dows 


An ideal way to make money with your hugerangeof | Fraxtorwindewe. Vovcte orto 
word-processor. Complete Marketing Pack musicalequipment _ | Voyetra $P Junior 
contains Clause Library, extensive Tutorial __ | Music Notation 


if i ms 
and Business development plans to produce ric cane hepcontice— Quick Score Deluxe 


professional quality Wills structured accord- Peck _ |Sound Cards 
ing to legal precedent. DOS and automated eat 

Windows versions - both supplied with fully 

bound example Will and other advertising 

materials. The whole system 


costs £75 - the price of apair = Other Software 
of Wills - the only computer Dp ‘Righotreet = Bent le-e-hox Au aranioing and 


F ; accompaniment. GS MIDI 
disc to pay for itself time X-OR 1.1 Universal editor/librarian 
and time again! 


eo Wave Windows based sound editor. 
Phone or write for Free Information Pack F ‘ ess Phone for best prices. If 
Legacy UK, (Dept D), PO Box 200 Tel: 0442-217624 | youdon —— you 
Shrewsbury SY11ZZ_ Tel: 0938 570428 a: want callenywey 
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Linux Base system £15 
Linux extras £19 
X-Windows (X11R35) £19 
Comprehensive System £45 
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The Public Domain & Shareware Library 
Winscombe House, Beacon Rd, 
Crowborough, East Sussex TN6 IUL 
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BLACK REFILL TWIN PACK 
BLACK REFILL TEN PACK 


COLOURED REFILL TWIN PACK 
(YELLOW, MAGENTA, CYAN, BURGUNDY, GOLD, GREEN) 


DESKJET CARTRIDGE (NEW) £13.99 
BJ10e CARTRIDGE (NEW) £18.78 
PRICES INCLUDE VAT AND P&P 


RECHARGED TONER 
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HP/CANON SERIES I, II, III 

EXTRA LIFE £35.00 + VAT 

IBM 401ie HIGH YIELD £70.00 + VAT 

HP III Si £60.00 + VAT 

HP/CANON SERIES I, I (BLUE,RED) £60.00 + VAT 
INCLUDES DELIVERY 


Linux is a System V like Unix OS for 
386 and 486 PCs. C compilers, 
programming libraries, the latest Unix 
user-interface all available. Fully 32- 
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networking. Prices exclude VAT. Call 
for our catalogue. Competitive 
hardware & PD software available 
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be bent round fixed points to form 
smooth curves while drawing yournext 
blueprint of a super space rocket or 
whatever. 

Realising that if a curve-generator is 
good enough for engineers, it’s prob- 
ably good enough for CAD work too, 
Bezier set about simulating the forces at 
work on a spline and created a short 
function which could be used to plot 
the bending ofa spline between points. 
The Bezier Spline (or B-Spline) algo- 
rithm has been gaining credibility ever 
since, and is now found on most self- 
respecting CAD systems. Just as in engi- 
neering, the spline is now a standard 
tool for most programmers. 

In fact, so prevalent is the Bezier 
curve in computing that it’s easy to 
forget just how it works, and why, some- 
times, it does not. The trouble is that 
while intuitively the Bezier curve is 
easy to use, the mathematics behind it 
can appear rather intimidating. Most 
programmers just use it to magically 
generate curves. 

Generally speaking, Beziers are cre- 
ated like this: two points are chosen by 
the user and a curve is drawn between 
them. The nature of this curve is set by 
two more control points which lie away 
from the curve but appear to act like 
magnets, drawing the bends ofthe curve 
towards them. By twiddling around 
with these, the programmer can create 
a suitably pleasing curve without wor- 
rying too much how it is generated. 

For a programmer curious to know 
exactly what the control points do to 
the curve, the formal definition ofthe B- 
Spline (given, as best it can be, by the 
formula in the SIMPBEZ function) is 
little help. The principle is solid enough: 
the line is drawn by varying t between 
0 and 1, and pulling out x and y values. 
(P1X,P1Y) and (P2X,P2Y) are the start 
and end points; (C1X,C1Y) and 


(C2X,C2Y) the control points. However, 
it is still not clear what is happening in 
this function, apart from it being a black 
box curve-making machine. 

Perhaps a better explanation is the 
diagram in Fig 1, for which I’m in- 
debted to Donald  Knuth’s 
METAFONTS. The illustration de- 
scribes a way of creating a B-Spline 
curve which relies on geometry rather 
than algebra. It’s not the most fun way 
to draw curves, but it does provide a 
more intuitive understanding of con- 
trol points. I think it’s worth doing just 
the once, if only to get the feel of how 
Bezier curves work. 

Start with the four points. Draw three 
lines: one between the first end point 
(P1) and the first control point (C1), one 
between C1 and C2, and one between 
C2 and P2. Now draw a line between 
the midpoint of C1C2 and the midpoint 
of P1C1, and between this point and the 
midpoint of C2P2. Now draw a line 
between the midpoint of thesetwo lines. 
In Fig 1, this gives a rough outline of the 
Bezier curve. If you need it more exact, 
the procedure can be repeated on the 
new set of points created. This has been 
done in Fig 1. 

Given this crude understanding of 
how the Bezier curve reacts to changes 
in its control points, you might like to 
rewrite Listing 1 to draw Bezier curves 
on your own machine. The language is 
AMOS, an Amiga-based Microsoft de- 
rivative (keeping true to an early Low 
Level promise not to forget the 68K 
machines). 

As it is the mathematics that are 
important, the program should be easy 
to convert. The ‘#’ sign after each vari- 
able name signifies floating point; most 
other versions of Basic default to this, 
so you may dispense with that symbol 
if you wish. The Bezier curve is imple- 
mented in two ways. One, SIMPBEZ, 
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. manipulates formulae, simplifies, 
integrates, differentiates, symbolically, 
solves equations, matrices, vectors etc. 


. E trig, exp, hyper, Bessel etc funcs 

. has 2D & 3D plotting, and 
recursive and iterative programming 

. is easy to use menus, split 


screen command and graphics. Now 
available for HP95LX Palm top. 


3-5 Cynthia Street, London N1-9JF 


CopyControl 


The ultimate Copy Protection System 
© User friendly & secure 


© Beats all bit copiers & disassemblers 

© Supports networks, backups, disk caches, windows etc. 
Remote changing of parameters 

@ No add-on hardware or special disks required 

© Control where, when & how often your programs are run 
© Compatible with all IBM MSDOS compatible computers 


For more information contact 


MICROCOSM LIMITED 


17 Cranbrook Road, Bristol BS6 7BL, UK 
Tel: (0272) 441230 Fax: (0272) 427295 


TIMELORD 


The simple-to-use Time and Materials 
recording package for Consultants, 
Solicitors, Architects etc; from the self 
employed to corporate multi-user........ 


Including DOS version, DOS package with PIF 
for Windows, TSR routine for on-line 
recording of phone calls, and a Load routine 
to pick up work from your portable, Timelord 
will improve your charge-out procedures and 
simplify your work recording, providing 
simple analysis of your profitability without 
hours of work. 
£95.00 Inc VAT 
Stubbs Systems 
1, Old Rectory Gardens,Godalming, 
Surrey GU7 1XB 
0483 419945 
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reflects the mathematical derivation 
of the curve and simply plugs a value 
of t into the Bezier function. CASBEZ, 
uses a faster recursive procedure based 
on the geometric derivation (an algo- 
rithm known as deCasteljau): Fig 1 
was drawn by a simple adaptation of 
this routine. 

Both are self-contained procedures 
that take the x and y coordinates of 
each of the four points as their argu- 
ments. The remaining code at the start 
of the program is a test routine which 
uses the mouse to move the four points 
around the screen, allowing a variety 
of Bezier curves to be drawn. A little 
tougher to translate line-for-line from 
Amigan, that; you are probably better 
advised to implement the spirit of this 
routine rather than get waylaid in the 
implementation’s idiosyncrasies. 
Non-adventurous PC types with e- 
mail facilities can contact me 
(dobrien@cix.compulink.co.uk) for a 
Visual Basic listing that achieves the 
same results. 

The most obvious thought to strike 
anyone playing with this short rou- 
tine is just how many smooth curves 
can be chosen that pass through the 
same two points. Most of the curves 
generated by the Bezier routine look 
fine. Take three points, however, and 
the rules become more stringent. It’s 
easy to join two Bezier curves together: 
simply share one end point between 
them. It’s a great deal harder to make 
the two curves join smoothly. Often 
(and again, experimentation is the best 
way to test this) a kink will form in the 
curve about the middle-point. 

The B-Spline routine is limited toa 
simple, two point curve. Moreover, 
there does not appear to be a simple 
way of extending it for more points. 
Mathematically, a logical route would 
be to extend the order of the function 
and make it a quadric. But this eats up 
computing resources, as it is rather 
shortsighted. What happens with a 
curve that has 99 points? Geometri- 
cally, it can be shown (don’t ask me to 
right now) that one Bezier curve can 
bend smoothly into the next curve, be 
continuous with it, ifC2 and P2 of the 
first curve and C1 of the second can be 
made to lie on the same line. 

That’s a useful rule when using 
Bezier curves manually to form shapes, 
but a little hard to implement auto- 
matically in a many-pointed curve 
drawing routine. Remember, we are 
attempting to find a function that will 
be able to draw a smooth function 
curve around a number of maxima 
and minima — up to five in the 
WALK.BAS program. The Bezier func- 
tion, while tried and tested, is not up 


Listing 1 
Screen Open 0,640,200,16,Hires 
Screen 0 


CAS_CUTOFF=3 : Rem deCasteljau Variables 


CAS_LEVEL=0 
Dim X(3),Y¥(3) 
MP=0 

Do 

Gr Writing 2 


For I=0 To 3 : Text 0,0,Chr$(48+MP) 


Next 
Repeat 
M$=Inkey$ 


If M$=”" “ Then Cls : For I=0 To 3 : Text X(I),¥(I),Chr$(484+I) : 


Next 
If M$<>”” and Val(M$)<4 Then MP=vVal (M$) 
Hard (X (MP) ) 


MX=X Mouse 

MY=Y Mouse 

MX=X Screen(0,MX) 

MY=Y Screen(0,MY) 

Text X(MP),Y(MP),Chr$(48+MP) 
Text MX,MY,Chr$(48+MP) 
X(MP)=MX : Y(MP)=MY 
MB=Mouse Click 

Until MB<>0 

Gr Writing 0 


: M$="3" 
: K Mouse=X 
: Y Mouse=¥ Hard(¥(MP) ) 


Rem SIMPBEZ[X(0),¥(0),X(1),¥(1),X(2),¥(2),X(3),¥(3)] 
Rem or CASBEZ[X(0),¥(0),X(1),¥(1),X(2),¥(2),X(3),¥(3)] 
Rem or BICUBIC[X(0),¥(0),X(1),¥(1),X(2),¥(2),X(3),¥(3)] 


Loop 
Rem 
Rem 


Procedure SIMPBEZ [AX#, AY#, BX#, BY#, CX#, CY#, DX#, DY#] 


Gr Locate AX#,AY# 
For T#=0 To 1 Step 0.01 


X#= (1-T#) * (1-TH#) * (1-TH) *AX#4+3*TH* (TH-1) * (TH- 


1) *BX#+3*TH*TH* (1- 


TH) *CX#+TH*TH#* TH*DX# 
Yi= (1-TH) * (1-TH) * (1-TH) *AYH#+3*TH* (TH-1) * (TH-1) *BYH4+3*TH*TH* (1- 
TH) *CY#+TH*TH*TH*DY# 


Draw To X#,Y# 
Next 
End Proc 


Procedure CASBEZ [AX#, AY#, BX#, BY#, CX#, CY#, DX#, DY#] 


Shared CAS_CUTOFF,CAS_LEVEL 
If (CAS_LEVEL) >=CAS_CUTOFF 


Polyline Ax#,AY# To BX#,BY# To CX#,CY# 


To Dx#, DY# 


Pop Proc 


End If 


Inc CAS_LEVEL 


AOX#=AX# : AOY#=AY# : DOX#=DX# : DOY#=DY# 
A1X#=A0X#+ (BX#-AOX#)/2 : A1Y#=A0Y#+(BY#-AOY#) /2 
B1X#=BX#+ (CX#-BX#)/2 : B1Y#=BY#+(CY#-BY#) /2 
C1X#=CX#+ (DX#-CX#) /2 : C1Y#=CY#+ (DY#-CY#) /2 
A2X#=A1X#+ (BLX#-A1X#)/2 : A2Y#=A1Y#+(B1Y#-A1Y#) /2 
B2X#=B1X#+ (C1X#-B1X#)/2 : B2Y#=B1Y#+(C1Y#-B1Y#) /2 
A3X#=A2X#+ (B2X#-A2X#) /2 : AZSY#=A2Y#+ (B2Y#-A2Y#) /2 
CASBEZ [AOX#, AOY#, ALX#, A1LY#, A2X#, A2Y#, A3X#, A3Y#] 
CASBEZ [A3X#, A3Y#, B2X#, B2Y#, C1X#, ClY#, DOX#, DOY#] 


Dec CAS_LEVEL 
End Proc 


to more than two points without a little 
tweaking of the control points. So how 
can we draw smooth, multiple-point, 
Bezier curves? 

It’s always a good idea, I’ve found, 
when stumped on a problem like this, 
to give up and try something else. In 
this case, the something else in ques- 
tion is not a million miles away. It’s 
another spline function, one I discov- 


ered quite by accident a year or so ago 
when digging around for sample Lisp 
code. One of the listings I uncovered 
was an implementation of this func- 
tion, the Bicubic Spline, perfect for our 
needs. Fundamentally, the function 
takes the same form. The formula is a 
function of t, where t varies from 0 to 1 
depending on how far down the curve 
you want to be, and x and y are pulled 
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386SX-25 2MbRAM/80Mb CALL 
38640 64KCache 4MbRAM/120Mb £CALL 
486SX25 64K Cache 4MbRAM/120Mb £CALL 
486-33 256K Cache 4MbRAM/120Mb £CALL 
486-50 256K Cache 4MbRAM/120Mb £CALL 


All Didar Systems include Small Footprint Desk Top Case, RAM 
and Hard Disk as Specified, 1.2Mb AND 1.44Mb Floppy Drives, 
SVGA Monitor, 1Mb SVGA Card, 2S/1P/1Games Port, Mouse, 
DOS 5.0, Windows 3.1, 102 UK Keyboard. 
One Year On-site Warranty. Call for other configurations. 
Call for latest prices 


SPECIAL OFFER COMPUTER PARTS 


OLIVETTI D33 LAPTOP 386sx33 
386dx 33MHz, 32K cache, 386sx33/64k 
60MbHD, 4Mb RAM, expand- THEY. 
ableto20Mb, builtin mouse 386dx40/64K......... 
pad, 32 greyscale VGA 486dx33/256k 
detachable screen, battery 486dx50/256k 
Trident 1Mb SVGA 


pack and charger PLUS 
FREE OLIVETTI HARD CARRY ; 
SVGA 14” Monitor ....£195.00 
102 Keyboard 


CASE WORTH OVER £120 
£999.00 
1Mx9SIMMS 


LASER PRINTERS 
OLIVETTIPG 404 = £549.00 
HP4* £1079.00 
MT904 Plus* £675.00 
RICOH LP 1200* £689.00 

OTHER PRINTERS 
CANON BJ10EX £185.00 

* ONE YEAR ON SITE WARRANTY * 

Call for prices on: Dot Matrix, Inkjet and other 

makes of printers, Motherboards, Drives, 


Memory, Monitors, Tape Streamers, 
Network Products, UPS and Cards. 


FREE DISK STORAGE BOX* 


3M (UNBRANDED) SONY/KAO/FUJI/DYSAN etc 


BOX OF 10 25 50 100 500 
3¥2"DS/DD 5.45 9.90 1885 35.65 156.25 
3¥2'DS/HD....9.25 19.80 36.45 69.70 303.60 
5¥%4"DS/IDD 4.15 645 11.50 21.85 101.20 
5¥s"DS/HD 7.40 12.10 21.40 39.75 178.25 


Call for pay discounts 
SPECIAL OFFER: 
DYSAN 312" DS/HD PREFORMATTED £9.00 PER BOX 
(Min 5 boxes) ALL DISKS 100% ERROR FREE 


5¥4" DISKS INCLUDES ENV/LABELS/WP 
31" DISKS INCLUDES LABELS 


DISK STORAGE BOXES 
De-luxe range: 
3%" 10Cap 40Cap 80Cap100Cap 54” 10Cap 50Cap 100 Cap 
1.50 3.95 4.45 4.95 1.60 4.10 4.60 


Stackable file drawers: 3%” 100 Cap £8.50. 54” 120 Cap £9.50 
A ah 


TONERS 
uP 


We also stock listing paper, fax 
rolls, copy holders, printer 
ribbons, cables & stands, 
mouse mats, mice, monitor 
stands, cleaning products, etc. 


PLEASE CALL FOR OUR COMPETITIVE PRICES 
“Deluxe disk storage box of your choice, absolutely free with every £100 
(net) worth of purchases or 10 boxes of branded disks. Offer only valid when 
this advert enclosed with payment and purchases made at advertised/ 
Catalogue | aoa E&OE. All prices exclude VAT @ 17.5% 
arriage from £3.20 + VAT per order. 

Next day delivery £9.00 + VAT at 17.5%. 

Hardware delivery POA. 


Cheques and PO's to: 


COMPUTHINGS EN 
63 Richings Way, Richings Park, Iver, 
Bucks SLO 9DB 
Tel: 0753630438 Fax:0753650009 
OPEN MON-SAT 8AM to 8PM _ ret:xxx9303 


SHAREWARE SETS 
FOR WINDOWS 


Buy one set for £12, two sets for £22, 
three or more sets for only £10 each 
ADD £1 POSTAGE PER ORDER (£2 EUROPE). 


Updated for 1993; even better value. 


EACH SET COMES ON SIX 3.5 INCH DISKS 
(1.44MB DISKS FOR THE GAMES + UTILITIES PACKS) 


WINDOWS GAMES PACK: OVER 160 GAMES! 


Contains over 160 games and leisure programs for 
Windows 3/3.1. Titles include Arachnid, Atmoids, 
Attaxx, Backgammon, Bang Bang, Chaos, Croquet, 
Checkers, Chomp, Fractal Paint, Fireworks, Fish, 
Football, Hearts, Islands, Jewel Thief, Klotz, Kye, 
Lunar Lander, Magic, Mines, Pool, Screenpeace, 
Snakes, Spacewalls, Sage, TilePuzzle, Windmaze, 
PipeDreams, Taipei, Warhead, Winfract, Winchess, 
Winjack, Winpoker, Wintrek .......... and many more! 


WINDOWS UTILITIES PACK: OVER 150 PROGRAMS 


Contains over 150 general programs including: 
Alfred, Aporia, Almanac, Active Life, Icon Draw, 
Icon Tamer, Icon View, Metz Utilities, Big Desk, 
Back Menu, Filefinder, Halfshell, Winzip, Winbatch, 
Organise, PaintshopPro, Trash Manager, Winspell, 
KwikDraw, Snagit, Winpost, Winedit, Dozens more 
utility programs, over 300 icons, fonts, bitmaps, 
picture file utilities. Over 10MB of programs for 
Windows 3/3.1 when unarchived. 


TRUE TYPE FONTS VOL 1: OVER 120 FONTS! 


Over 120 different shareware and public domain typefaces 


AlbatrossArchitect Arctic BfackChanceryCaraway 
CrtWrighCaruP epreaCrillee<wie7ontQaytona 
Dubiel LINTFONTF Tank TimesGraphicl ightHfeideberg 
\ibbywcuiptlowarEastlLowerWestSidel unpasky 
MachumaineO sw aldG rey PixiePoloSemiRivrside 


RudelsbergSeript Jantinformal cy CeJoulouselgutree 
Upped st dU tra Black (AGERE and morel! 
All are in the new Windows 3.1 TrueType format!! 
TRUE TYPE FONTS VOL 2: OVER 120 MORE! 
A full printout is included with each collection 
PrdromedabLACKADDERBRIGH TONGairoCamber lic 
ChopinDebussyDiego UAW OODEnglishBown 
FaktosFlemishGismondabkecwWar'SK ASTELLAR 
Koshgariank OOP. ilithlighf WARBOLO MAZAMARE 
MiddictonMQULTIFORM Paganini PCMiraPianost 


Rent ofSinoteSchwarzwald/ ANANCAL /THGMONITORS 
Noo X ov reSeo ce A¥ide0 VreoFot ZALESKI and morel! 


Easily installed; complete instructions supplied 
a a eee eee 


I Please send me the following Windows sets 4 
I Q) WINDOWS GAMES Q) WINDOWS UTILITIES | 
| Q TRUETYPE FONTS 1 Q TRUETYPE FONTS 2 | 
I ) l enclose a cheque/postal order for £ ...... 
I Q Debit my credit card number 


QO00 QO00 OO00 OO00 
| Expires QU OU Signature: 


1 LL eeeteeeranaistaiipc B58 ORE AEE TERRE RTO Oc TELS 


I PGUPOSE | acti veassinisecsiepsesanascUninsnddcudeshedverntsaxes 


CS9302 


SHAREWARE CONNECTIONS 


wee ey TEL: 041-775:1937 en 
== LE) 49 BOGHEAD ROAD © 


KIRKINTILLOCH GLASGOW G66 4ED 


from MKS 


provides an easy, step-by-step introduc- 
tion to the operating system, file system 
and UNIX commands- 

¢ work with the Korn Shell 


« use UNIX shell scripts on your DOS PC 
* includes 160 UNIX commands 

* comes with tutorial software 

* requires 286 processor or better 

* complete package only £119.00 plus VAT 


from your/Mf KS specialist 
System Science 
3-5 Cynthia St, London N1-9JF 
(071) 833-1022 


EASYSIGN 


An easy way to create signs, posters and 
banners on any |BM-compatible PC. Letters 
print at any size, restricted only by printer 
& page sizes. Prints on Epson-compatible 
dot-matrix printers. 

Available on 32" or 5%" disks. 


At only £19.95, how can you 


afford NOT to have it? 
(price includes VAT and postage) 


also, powerful software for cutting 
self-adhesive vinyl. Ring for details. 


WIGHT SCIENTIFIC 


9 Blackheath Hill, Greenwich, London SE10 8PB 
Tel: 081-694 9293 Fax: 081-694 8946 
Mastercard & Visa accepted. 


SEND YOUR COMPUTER oo 
TO MEDICAL SCHOOL! (a= 


12,000 copies sold worldwide 


Sending your computer to medical school 
helps you stay healthy 
“Doctor Schueler's Home Medical Advisor” gives you an enor- 


mous amount of accurate medical information directly from your 

PC in a few keystrokes. It is a complete home medical system 

developed by United States Board Certified Physicians. It com- 

prises a set of integrated files giving a complete medical encylo- 

pedia with full colour graphics and straightforward explanations 

of medical terms. 

@ Fully integrated and easy to use, 

@ Print or export information onto Hard Disk. 

@ Uses a surgery simulation with over 750 colour graphics to 
diagnose over 500 common ailments and diseases. 

@ Allows you to identify and see the side effects of over 1000 
prescriptions and over the counter drugs. 

@ See how to treat over 150 injuries. 

Also describes standard medical tests, gives treatment for house- 

hold ingestions, and health and diet advice to keep illness at bay. 


ONLY £59.95 inclusive 
PC compat 512K RAM HD EGA/VGA 


Distributed by — AKE Technology, 9 Peel Avenue, 
Retford, Notts, DN22 7JU. Cheques/PO only. 
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Listing 2 
Rem FOR BICUBIC, Insert: 
Rem X(0)=X(1) : ¥(0)=¥(1) 
Rem X(1)=X(2) : ¥(1)=¥(2) 
Rem X(2)=X(3) : ¥(2)=¥(3) 


ae HANDS ON - LOW LEVEL 


rem just before the Loop command in Listing 1 
Procedure BICUBIC[I0O,J0,I1,31,12,32,13,33] 


WEIGHT#=0.5 
PIXRES#=2.5 
Gr Locate I1,J1 


DIST#=Sar ( (11-12) * (I1-1I2) +(d1-32) *(J1-g2) ) 


If (PIXRES#<DIST#) 
75=J1-J0 
I5=I1-10 
J3=J32-33 
I3=I2-13 
DT#=PIXRES#/DIST# 
COUNT=Int (1.0/DT#) 
While (COUNT>1) 
Dec COUNT 
T1#=DT#* COUNT 
T2#=1-T1# 
TS#=WEIGHT#*T1#*T1#*T2# 
T3#=WEIGHT#*T1#*T2#*T2# 


J=I1*T1H4+I02*T2#4+T3*T3 #475 * TSH 
T=I1L*T1#+12*T2#+13*T3#+15*T5# 


Draw To I,J 
Wend 
End If 
Draw To I2,J2 
End Proc 


out from the equation which, for the 
curious, is 

x = xl*t + x2*(1-t) +W*(x2- 
cx2)*t*(1-t) *(1-t) +W* (x1- 
cx1) *t*t*(1-t) 
for the x coordinate, where x1 and x2 
are the two end point x coordinates 
and cx1 and cx2 are the two control 
point x co-ordinates. It’s the same for- 
mula for y, of course. 

Geometrically? Well, I’m afraid I’m 
no expert, but the illustration in Fig 2 
should provide the same rule ofthumb 
as Fig 1, if not the demonstrative rig- 
our. Roughly speaking, the curve gen- 
erated between P1 and P2 tends 
towards C1P1 at P1 and C2P2 at P2. 
The strength of the tendency at P11 is 
determined by the distance of C1 from 
P1. The further away the centre point 
is, the closer follows the vector de- 
scribed by the line. 

Inother words, the curve P1P2 tends 
to act as a smooth extension of the 
lines C1P1 and C2P2. Given this, it 
doesn’t take long to realise that a 
smooth, multiple-pointed curve can 
be drawn by making the control points 
of the present curve point to the end 
point of the previous curve and the 
beginning point of the next. To draw a 
smooth curve made up of the points 
P(1)...P(n), starting with acurve drawn’ 

y a procedure call 
Bicubic(O,P(1),P(2),P(3)), followed by 
aBicubic(P(1),P(2),P(3),P(4)) and end- 
ing with Bicubic(P(n-2),P(n-1),P(n),O). 

Listing 2 gives an implementation 


of the Bicubic Spline algorithm, again 
for AMOS. By entering in the few lines 
before the main procedural definition 
after the DoSpline call in Listing 1, the 
code can use the test routine to draw an 
infinitely long smooth spline. 

There are a few problems with the 
Bicubic Spline routine. Firstly, speed. 
While it matches a simple mathemati- 
cal implementation ofthe B-Spline, itis 
noticeably slower than the deCasteljau 
Bezier curve. Secondly, it is not always 
easy to know what values to put at the 
very ends of a long Bicubic spline (the 
values given as ‘O’ in the above para- 
graph). The obvious solution is to sim- 
ply repeat P(1) at the beginning of the 
curve and end with P(n). Other, more 
complex algorithms — notably, a rou- 
tine known as the de Boor algorithm — 
provide smoother end conditions 
known as not-a-knot end points. 

If the spacing between points is too 
wide, the function has a tendency to 
overshoot slightly around P1 because 
the degree of closeness between the 
beginning of the P1P2 curve is propor- 
tional to the distance between P1 and 
C1. This can be tweaked by changing 
the WEIGHT# variable in our algorithm. 
Actually, this aspect of the Bicubic 
Spline is both a bug and a feature. In 
curves where a large distance implies 
great speed (for instance, in the regu- 
larly-sampled position of an object fly- 
ing through the air), the overshooting 
reflects the natural curve followed by a 
fast-moving but decelerating object. The 
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Bicubic is best suited to closed out- 
lines and frequently sampled curves. 

Having wandered far from the origi- 
nal problem, it might be a good idea to 
recap here. Curves are most often im- 
plemented on computers in order to 
draw appealing graphical representa- 
tions of natural objects, although they 
can also be used in curve-fitting — the 
attempt to interpolate between a set of 
sampled points using a smooth, natu- 
ral function. I haven’t had the space to 
cover this aspect of curve simulation 
but I hope it should be clear how the 
Bicubic Spline routine can be adapted 
to fix the clumsy nature of the present 
WALK.BAS interpolation routine. 

When it comes to drawing a short 
(two point) curve ona graphics screen, 
Bezier curves are fast and aesthetically 
pleasing. However, for larger and more 
complex curves, the Bicubic Spline 
routine is a useful general purpose 
routine. And while it has problems 
with speed, end points and a mild 
tendency to overshoot, it still beats 
Bezier into a cocked hat. 

Two books are worth investigating: 
the alliterative Bartels, Beatty and 
Barsky’s An Introduction to the Use Of 
Splines in Computer Graphics (Menlo 
Park, California. Morgan Kaufman, 
1986). That old bible Numerical Reci- 
pes in C (Press, Flunney et al. Camb. 
Univ. Press, 1988) covers curve-fitting 
reasonably well, too. The algorithm 
given here for Bicubic Splines is taken 
from Lauzzana and Penrose’s Imple- 
menting Bicubic Splines in Allegro Lisp 
(Dr Dobb’s Journal, August 1990). 


QBasic users retaliate 
Finally, a large apology to the irate 
QBasic users who wrote and called in 
about December’s piece on that lan- 
guage’s shortcomings. The article 
claimed, at various points, that the 
QBasic supplied with DOS 5.0 could 
only cope with 160K of program and 
data, that it did not have local vari- 
ables, and that it lacked an ON...EVENT 
handler. None of these statements is 
true. QBasic can cope with 64K of 
program persubroutine, automatically 
assumes local variables, and has a full 
event handler (although as Neil Lomas 
of Crewe noted, it does seem to have 
problems handling timed events dur- 
ing INPUT statements). If it is any 
compensation, the main problem 
pointed out by the piece (the lack of a 
CALLINTERRUPT command) and the 
solution provided are entirely valid. 
If you have any suggestions, correc- 
tions or additions to Low Level, don’t 
hesitate to write (via the PCW offices) 
or email me _ direct on 
dobrien@cix.compulink.co.uk. 


